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Abstract — This paper studies joint beamforming and power 
control in a coordinated multicell downlink system that serves 
multiple users per cell to maximize the minimum weighted 
signal-to-interference-plus-noise ratio. The optimal solution and 
distributed algorithm with geometrically fast convergence rate 
are derived by employing the nonlinear Perron-Frobenius theory 
and the multicell network duality. The iterative algorithm, 
though operating in a distributed manner, still requires instan- 
taneous power update within the coordinated cluster through 
the backhaul. The backhaul information exchange and message 
passing may become prohibitive with increasing number of 
transmit antennas and increasing number of users. In order to 
derive asymptotically optimal solution, random matrix theory is 
leveraged to design a distributed algorithm that only requires 
statistical information. The advantage of our approach is that 
there is no instantaneous power update through backhaul. More- 
over, by using nonlinear Perron-Frobenius theory and random 
matrix theory, an effective primal network and an effective 
dual network are proposed to characterize and interpret the 
asymptotic solution. 

Index Terms — Power control, coordinated beamforming, max- 
min duality, effective network, large system analysis, multicell 
network, nonlinear Perron-Frobenius theory, random matrix 
theory. 



I. Introduction 

TO benefit from the available and increasing spatial de- 
grees of freedom, multicell networks exploit different 
forms of intercell cooperation to operate the system in an 
interference-aware manner JT], (Q- Due to practical constraints 
such as limited feedback |3|-|5| and the finite capacity of 
the backhaul |6)-|8), beamforming level coordination and 
efficient power control strategies are favored over data level 
cooperation and nonlinear precoding approaches [9], [ 1 1 
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to effectively scale up the system performance. Considering 
these practical constraints, two characteristics are appealing 
to joint beamforming and power control algorithms design: 
distributed computation and fast-convergent algorithms with 
low complexity. The desired distributed feature addresses 
system scalability, and the distributed algorithm only relies on 
local channel state information (CSI) which can be obtained by 
uplink measurement in a time division duplex (TDD) system 
or through user feedback in a frequency division duplex (FDD) 
system. On the other hand, simple algorithms possessing fast 
convergence rate are attractive in that they reduce the message 
passing overhead and alleviate the finite backhaul constraint. 

The algorithm design is intimately related to the system 
performance metric of interest. Different system performance 
metrics reflect different design priorities. One common ap- 
proach is to maximize the sum rate of the system. How- 
ever, due to the non-convexity of the problem, numerically 
finding the optimal solution is challenging and the design of 
distributed algorithms that can compute the global optimal 
solution efficiently is still open, e.g., see (TJ, lfTTI -lfT9l and the 
references therein. It is known that two specific problem for- 
mulations admit global optimal solutions: the transmit power 
minimization subject to signal-to-interference-plus-noise ratio 
(SINR) constraints, and the maximization of minimum SINR 
subject to power constraints. The former problem whose 
priority is energy saving has been addressed extensively in the 
literature and efficient algorithms have been proposed for both 
the single cell and multicell systems [ 20 1— [ 3 1 1. The analysis of 
single cell downlink relies on the well-known uplink-downlink 
duality |23l . [25|-[28| which is readily interpreted by the 
Lagrange duality in convex optimization. In |29l , the duality 
is observed for the MIMO multiuser ad hoc network setting, 
and in [30], the duality is extended to the multicell setting. 

The literature for the latter problem which aims to enforce 
the fairness level of the system is comparatively less. The 
max-min SINR problem was first addressed in 11321 using 
an extended coupling matrix approach, and a centralized 
algorithm was proposed in 11271 . which involves an increased 
dimension matrix computation. A reformulation of the max- 
min problem is analyzed in [331 by conic programming and 
a heuristic algorithm is provided. In 11341 . similar to [33], 
the max-min problem is tackled from the transmit power 
minimization perspective and a hierarchical iterative algorithm 
is proposed. Recently, the problem was studied in lfl5l using 
a nonlinear Perron-Frobenius theory H351 . and a distributed 
algorithm was proposed that exhibits the distributed power 
control (DPC) structure in [20 1 . The DPC-like structure is 
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independent of parameter configuration, thus enabling the 
application of the power control module in [20 1 already used 
in practical cellular systems. The approach |fj"5l is extended 
to the MIMO downlink in f36l wherein the convergence of a 
heuristic algorithm in 11331 is proved, and a power optimization 
problem under multiple power constraint is analyzed in |[37l . 
The optimization of the egalitarian fairness i.e., max-min, 
performance metric also has intimate relationship with other 
important wireless network performance metric optimization 
problem, e.g., the weighted sum rate maximization [38|-[40|. 
The sum rate maximization is nonconvex and NP-hard, and 
the implication is that fast egalitarian fairness algorithms can 
be leveraged to solve this nonconvex problem with global 
optimality guaranteed under special cases. Herein, we firstly 
extend the analysis in |[T5l . [36] to the multicell setting with 
multiple serving users per cell. The duality between primal 
and dual network is derived and characterized by the Perron- 
Frobenius theory. A distributed algorithm is also proposed 
which possesses geometrically fast convergence rate. 

The designed algorithm, though converging to the optimal 
solution, requires instantaneous power update within the coor- 
dinated cluster through backhaul. This instantaneous informa- 
tion exchange may become prohibitive when the number of 
transmit antennas at base station as well as the serving users 
per cell grow large. In such emerging large-scale multiple 
antenna systems |RD - Il44l . the backhaul capability may turn 
into the bottleneck. In order to alleviate this problem and 
to enable simplified design that utilizes only the statistical 
channel information, additional tools from random matrix 
theory 11451 . (46 1 are to be leveraged. The large system 
analysis for linear receiver design in the uplink was initiated 
in l47ll . and the notion of effective interference and effective 
bandwidth was proposed. In ll48l . asymptotic analysis for 
the transmit power minimization problem is carried out. The 
approaches in [49 1 and [50] decouple beamforming and power 
control by assuming zero-forcing or regularized zero-forcing 
beamformers 11511 . The analysis in ll52l examines the max- 
min SINR problem from the transmit power minimization 
perspective, and compares several cooperation strategies by 
assuming a two-cell model with homogeneous channel setting. 
The analysis is extended in [53 1 for the minimization of the 
maximum power problem with homogeneous channel setting. 
In this paper, we perform large system analysis for the max- 
min SINR problem in a general multicell setting. Utilizing 
tools developed from random matrix theory, the deterministic 
equivalents ll46l . lf54l for the dual network SINR and for the 
primal network SINR are established. These asymptotic ap- 
proximations are used to compute the asymptotic power which 
only relies on statistical channel information. Intuitively, in a 
large-scale multiple antenna system, the optimal powers for 
different users would approach different deterministic values 
and the obtained power can be utilized for optimal beamformer 
design with local CSI. Moreover, by using nonlinear Perron- 
Frobenius theory and random matrix theory, we observe an 
effective network for the dual network and an effective net- 
work for the primal network, which capture the characteristic 
of the power control effect in the large system setting. The 
established effective network is further leveraged to provide a 
distributed algorithm with fast convergence rate. 



To summarize, the contributions of this paper are three- 
fold: 1) analysis and algorithm design for joint optimal beam- 
forming and power control in a finite multicell system to 
maximize the minimum weighted SINR, 2) the established 
effective network to characterize the algebraic structure of 
the power control problem in the large system setting, and 
3) low complexity algorithm design which requires no in- 
stantaneous backhaul exchange. All these contributions lead 
to efficient methodologies to design algorithms for the large- 
scale coordinated multicell downlink. The paper is organized 
as follows. Section [n] presents the system model. The finite 
system analysis is provided in Section [Till Section [TV] carries 
out large system analysis and derives the asymptotic solution. 
Numerical results are presented in Section |V] Finally, Section 
[VI] concludes the paper. 

Notations in this paper are presented as follows. Boldface 
upper-case letters denote matrices, boldface lower-case let- 
ters denote vectors, and italics denote scalars. The Perron- 
Frobenius eigenvalue of a nonnegative matrix F is denoted as 
p(F). Let x(F) and y(F) denote the Perron (right) and left 
eigenvectors of F associated with p(F) respectively. Tr(A) 
denotes the trace of the matrix A, and diag(a) denotes the 
diagonal matrix having the vector a on its diagonal. Let 
(/(a)) m denote the mth element of a function vector /(a). 
Let aob = (ai&i, • • • , om&m) T (the Schur product). Let C, 
R+, and R++ represent the set of complex numbers, the set of 
nonnegative real numbers, and the set of positive real numbers 
respectively. Let (-) T and (•)T denote the transpose operation 
and conjugate transpose operation respectively. || • || denotes 
the Euclidean norm for vectors and spectral norm for matrices, 
and — ^ denotes almost sure convergence. 

II. System Model 

Consider a coordinated multicell downlink formulated by J 
coordinating base stations utilizing the same carrier frequency. 
Each base station is equipped with N transmit antennas and 
serves K users simultaneously. Herein, the focus is on the 
base station side interference coordination, and each user is 
assumed to have a single antenna. The received signal yj % 
for user k in cell j is written as 

j 

where h/ j fc £ C*** 1 denotes the channel vector from cell 
I towards user k in cell j, Xj <E C Nxl is the transmitted 
signal vector of cell I, and z^k characterizes the additive white 
noise effect and any intercell interference not included in the 
coordinated cluster for user k in cell j, which is distributed 
as CAA(0, <7j, fe ) with cr hk G 

Linear beamforming strategy is assumed at the base station, 
and thus the transmit signal vector Xj for cell j can be 

expressed as Xj = Y,k=i x j,*> = J2k=i y^ s j,fc u .7,fc> where 
Xj,k S C Nxl represents the signal intended for stream k of 
cell j, Sj t k and denote the information signal and the 
transmit power for that stream, and Uj j. € C Wxl denotes 
the normalized transmit beamformer for user k in cell j, i.e., 
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l u j,fc|| 2 = 1. The SINR for user k in cell j can be written as 



r™A SINR PN = 



^rl h L,fc u j.fel 



(2) 



where the superscript (-) PN represents the primal downlink 
network. Let uijk denote the weight associated with p,jj~ for 
user k in cell j illustrating different power prices, and denote 
(3j } k as the priority factor associated with r PN , for user k in 
cell j demonstrating diverse service priorities. Then the max- 
min problem under weighted sum power constraintQ can be 
written as follows 



maximize 



ft,* 



subject to Vir,. A ^- /'. p i)fc >0, \\u^ k \\ 2 = l 



variables : Pjk, Uj k ■ 

(3) 

The problem (0 appears non-convex at first, but can be 
transformed into a second-order cone program tl55l by ap- 
plying methods similar to that in [33 1, which admits a 
global optimal solution. However, employing standard convex 
optimization methods to find the optimal solution typically 
requires centralized computation and incurs a fair amount of 
parameter tuning and message passing overhead that may not 
be practical in wireless networks. Thus in Section Hill we will 
employ nonlinear Perron-Frobenius theory to propose DPC- 
like algorithm [20 1 that does not require parameter tuning 
and has geometrically-fast convergence rate. Then in Section 
II VI algorithms that are even simpler and more practical for 
systems with a large number of transmit antennas and users0 
will be presented by performing an asymptotic analysis. 

III. Finite System Analysis 

This section is devoted to finite system analysis when N 
and K are not asymptotically large. Section IIII-AI reformu- 
lates problem (0 to exploit its analytic structure. Section 
IIII-BI establishes the network duality via a Perron-Frobenius 
characterization, and provides a geometrically-fast convergent 
algorithm to compute the optimal solution. 

A. Problem Reformulation 

The problem formulation in ((3) essentially regards an in- 
terference network with JK users. However, the formulation 
in terms of the channel h^ fc and the link gain |h/ fc u i,i | 2 
does not easily lead to amenable analysis. In order to construct 
the JK x JK cross channel interference matrix, consider the 
matrix G e R'^ xJK with subscripts m and n, whose entry 
can be written as 



G m , n - |hf„ ut1 !m _Aifj u rtlrfLf jf 



(4) 



where |~-] and [ J denote the ceil and floor operation re- 
spectively. Thus the channel hj can be represented with 

'The weighted sum power constraint implies that multiple base stations 
form a coordinated cluster to jointly perform power control in order to achieve 
the desired fairness level for users in the cluster. 

2 The large system analysis with algorithm design for a single cell downlink 
is studied in |56|. 



subscripts m and n: h mj „ = h^n-^^m^ m _ K ^mj. Moreover, 
define the power vector p 6 M^ xl as p m = p|-m.^ TO _^-|_s>j, 
and the beamforming matrix as U = (ui,-- - , ujk) with 
u m = urm-| )m _j^|mj . The general formulation in (0 can 
be easily interpreted through two special cases: a) J = 1, K 
arbitrary and b) K = 1, J arbitrary. The former case refers to a 
single cell downlink with K interfering users, while the latter 
case corresponds to an ad hoc interference network setting 
with J transmitter-receiver pairs or a multicell setting with 
one user served per cell. By the formulation of G, the cross 
channel interference matrix, denoted by F £ M.+ Kx JK can be 
obtained by 



if m = n 
m , n , if m ^ n. 



(5) 



Similarly, the weight vector w e R^ xl , the priority vector 
(3 6 Ril xl , and the noise vector <x e ™ JKxl 



-+ 



can be 



defined by: w n 



w 



and 



c m = (7r»| iffl _K|«i. From the aforementioned mapping, if 
we denote the SINR vector as T PN e M^ xl with T,™ = 

pPN „„j .u- . r- m>JKxl 



and the auxiliary vector g E 



with 



A 

g = 



1 

Gi.i 



' G JK..JK f 

(0 can be readily reformulated as follows: 



then the optimization problem 



maximize 

subject to 
variables : 



J- W T 

N 

p,U. 



r™( P ,u) 
P<P, 



~ (diag(/3og)(iFp+ C r)) 

p>0, ||u m || 2 = l 



(6) 



It can be shown that solving (|6]l is equivalent to solving (fj). 
The compact formulation in (O introduces a nonnegative ma- 
trix diag(/3og) (F + (l/P)erw T ), whose algebraic structure 
helps in establishing the network duality and is pursued next. 

B. Network Duality and Algorithm Design 

The analytic structure in © is similar to the formulation in 
[36 1 for the single cell multiuser downlink scenario. In [36 1, 
the uplink-downlink duality is proved by a geometric pro- 
gramming formulation and the Lagrange duality. Herein, we 
provide a network duality interpretation for the max-min based 
multicell scenario via Perron-Frobenius characterization. 

For any given beamforming matrix U, a simpler optimiza- 
tion problem for © can be formulated by only optimizing the 
power solution. It is known that at optimality, the weighted 
SINR for different users are the same, and the weighted power 
constraint becomes tight 03]. Now if we explicitly make the 
dependence on U and denote the optimal weighted SINR as 
r*(U), then the optimal power solution satisfies [15|, ll36l : 



1 P*(U) 
r*(U) N 

(7) 

From (0, it can be shown from nonnegative matrix 



= diag(/3og(U)) (F(U) + (l/P)<rw T ) £M. 



theory [57 1 that 



N 



is the Perron (right) eigenvec- 



tor (up to a scaling factor) of the nonnegative matrix 



namely, p — 

^(^,ffO)(FCU) + (i/^) ) and r , (u) is related to ks 

" )) 



diag(/3 o g(U)) (F(U) + (1/P)c 



diag(/3og(U)) F(U) + (l/P)c 
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TABLE I 

Algorithm A: Max-min weighted SINR for multicell downlink 



. Initialize arbitrary p[0] G , q[0] G R++ Xl and u m [0] G 

C JVxl for m = 1.....JK such that ||u m [0l|| = l.Vm, 



iw T p[0] < P, and i<r T q[0] < P 
1) Update dual network power q[re + 1] 



9m [« + 1] 
2) Normalize q[/c + 1] : 



/3„ 



rDN( qM ,UM) 



NP 



g TO [«] Vm. 



q[/t + 1] <- q[* + 1]. 

er 1 q[re + 1] 

3) Update transmit beamfonning matrix U[« + 1]: 



[k + 1] 



u m = r^Ti 1 Vm. 



4) Update primal network power p [« + 1 

Pm[K+11 = (™l 

5) Normalize p[/e + 1]: 



U[« + l]) 



Pm[«] Vm. 



p[« + 1] 



/VP 



w T p[re + 1] 



p[«+l]. 



spectral radius by the following: 

T*(U) 



1 



p (diag(/3 o g(U)) (F(U) + (l/P)<rw T )) 



(8) 



Now, in order to establish the network duality, consider the 
hypothesized dual uplink network by construction. Denote the 
dual network transmit power solution q € R^* 1 as q m = 
<7rm-| im _#L»Bj, where with subscripts j and k denotes 
the reciprocal uplink transmit power for user k in cell j. Let 
the weight vector w in the primal network be the noise vector 
in the dual network, and conversely let the noise vector er in 
the primal network be the weight vector in the dual network. 
Then the max-min problem for the dual network given receive 
beamforming matrix U can be formulated as 



maximize mm 



r™( q ,u) 

0m 



(diag(j3og(U)) ( £ FT (U)q(U) +w) ) v 



subject to i(T T q(U) < P, q(U) > 
variables : q(U) 

(9) 

where the superscript ( ) DN denotes the dual uplink network. 
By leveraging the following properties of nonnegative matrices 
11571 : p(A) = p{A T ) and p(AB) = p(BA), the optimal 
solution for (O equals —r- — — — , — tvt- Com- 

H - J 4 p(diag(/3og(U))(FT(U) + (l/P)w -T)) 

paring with the optimal solution for the primal network in 
(JHJ, the network duality is observed by employing F T as 
the cross channel interference matrix for the dual network 
and reversing the role of w and er. Thus the optimal power 
solution given U is the left eigenvector of the nonnega- 
tive matrix diag(/3 o g(U)) (F(U) + (l/P)erw T ), namely, 

q*(U) _ NPy(a iae (f3o s (V))(F(V) + (l/P)™ J )) N , . 

~ CTT y(dia g (/3og(U))(F(U) + (l/P) CT wT)) " N ° te ttlat SmCe 

the network duality holds for any given U, the achievable SINR 
regions of the max-min problem are the same for both the 
primal network and the dual network. 



The motivation for establishing the dual network is to 
exploit the decoupled property of the receive beamformer 
optimization and to utilize the optimized received beamformer 
as the optimal transmit beamformer for each user. The optimal 
beamforming matrix U* depends on the power vector q, and 
for any given q, the optimal beamformer u^(q) can be 
obtained by 



<i(q) = argmin- 



%jtm N m , n m,n 



l rn LL m,m LL m,m 



(10) 

which can be readily solved and is known to be the minimum 
variance distortionless response (MVDR) beamformer which 
is given by: 

u m(q) 



t 

m.n 



l I)- 1 h, 



(11) 



Therefore, the optimal solution for the beamformers, the power 
of the dual network, and the power of the primal network 
can be written as: = u^(q*), q* = q*(U*), and p* = 
p*(U*). The optimal solution is of analytical interest. In order 
to derive a fast algorithm to compute the optimal solution in 
a distributed manner, we employ nonlinear Perron-Frobenius 
theory and propose the algorithm given in Table Q] referred 
to as Algorithm A for the multicell scenario. It exhibits the 
DPC-like structure as in lfT31 . ||36 I for the single cell scenario. 
The convergence property of Algorithm A is discussed in the 
following theorem. 

Theorem 1. Starting from any initial point q[0], p[0], and 
U[0], the q[n], p[«], and V[k] in Algorithm A converges 
geometrically fast to the optimal solution q*, p*, and U*. 

Proof: The proof is given in Appendix A. ■ 
Remark: Distributed algorithms utilizing only local CSI 
and requiring limited backhaul exchange are important for 
practical implementation issues. Algorithm A is distributed in 
the sense that the iterative update (step 1, 3, 4 of Algorithm 
A) can be independently performed for each individual user at 
each base station. In addition, each base station only employs 
local CSI, which can be directly obtained in a TDD system 
or acquired by user feedback in a FDD system. The normal- 
ization procedure (step 2 and 5 of of Algorithm A), however, 
requires a central computation of w t p[k] and er T q[K]. This 
procedure can be made distributed by gossip algorithms [58] 
and power update through the backhaul. 

Hitherto, an algorithm for computing the optimal solution to 
(0 is established. In Section |VJ we provide numerical results 
that support and confirm its fast convergence property. Further- 
more, with minimal parameter exchange and configuration, 
this algorithm is practical in a finite system. However, in a 
large-scale system when both and K become large, the in- 
stantaneous power update across the coordinated cluster limits 
its practical implementation. Therefore, a lower complexity 
algorithm is needed in large-scale systems and is studied next 
in Section |TV] 

IV. Large System Analysis 

This section is devoted to a large system analysis when 
both the number of transmit antennas N and the number of 
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serving users per cell K go to infinity while the ratio (load 
factor) lim ^ remains bounded, i.e., the notation N — >• oo 
denotes that both N and K become large, while lim inf |r > 
and lim sup ^ < oo. In this large-scale system setting, for a 
given channel realization, the amount of instantaneous power 
update through the limited backhaul can be impractically large 
and thus impacts the system performance. A key question 
is whether it is possible to design a non-iterative algorithm 
to compute the beamformer and still achieve some form of 
optimal egalitarian fairness. Herein, the optimality is in the 
asymptotic sense. This means that, if the power p and q 
in the large system converge to some deterministic values 
that only rely on statistical channel information, then these 
deterministic values can be a priori calculated, stored, and 
updated only when the channel statistics changq^ Thereafter, 
the beamforming matrix ought to be non-iteratively computed 
using these slowly updated power values and the available 
instantaneous local CSI. 

This idea of practical implementation for large systems 
will be studied by addressing two problems related to (0. 
Firstly, different users in the multicell network have potentially 
different weights, different priorities, different noise powers, 
and more importantly, different large-scale channel effects 
which may consist of path loss, shadowing, and antenna 
gain. Thus, to maintain the max-min fairness across users, 
the powers for different users would converge to different 
deterministic values in the large system setting. One key 
issue is to establish the asymptotic optimality for both the 
dual network power q and primal network power p. Another 
key issue is to design distributed algorithm to compute these 
deterministic values. 

In Section [TTT1 no specific channel models are assumed. 
Now for amenable analysis, the transformed notation using 
subscripts m and n will be still employed and the following 
channel model is further assumed: 



(12) 



where d m , n represents the large-scale channel effect and il- 
lustrates the statistical channel information. The h m n denotes 
the normalized CSI whose elements are independent and iden- 
tically distributed as CAf (0,1). This assumption corresponds 
to the practical setting where the antenna elements equipped 
at each base station are placed sufficiently apart. Herein, 
independent channel assumption is employed and the analysis 
with the general correlated channel model [71 1-[75 | is left for 
future work. Employing this channel model, the asymptotic 
analysis for the dual network and primal network is carried 
out in Section IIV-AI and Section IIV-B1 respectively. 

A. Asymptotic Analysis for the Dual Network 

The large system analysis for the dual network is examined 
first to derive the asymptotic dual network power, which is 

3 This idea relates to algorithmic developments in the context of and 
in support of the design of situational aware wireless networks |59|. The 
envisioned situational aware wireless networks adapt system parameters and 
algorithms design to the channel attributes (i.e., different types of channel 
information, various channel statistics reflected in different dimensions), user 
attributes (i.e., different user densities, fairness requirements, user mobilities), 
and system attributes (i.e., backhaul capabilities, large system or sparse system 
structures, energy efficiencies), which constitute the wireless environment and 
network situations, see [601—170] for examples driving this trend. 



utilized for beamformer design. One key step is to study 
the asymptotic behavior of the dual network SINR, whose 
expression is given by using the optimal MVDR beamformer 
as follows: 



r° N (q) = %hm,m \ ^ ^hm,nht,,» +W m I ] h 



N 



Q m dm . m {" f / V ^ Qn dm ,n f ft t I f w 

h m,m I fj h ™," h m,n + ™™ 1 h ™.m Vm - 

\ n^m J 

(13) 



Since each instantaneous CSI is random, the instantaneous 
SI N R in ( fl3l l is a random variable in quadratic form. Moreover, 
since the dual network power and large scale channel effects 
are diverse across users, if we define the random matrix H m as 
H m Hj„ = Y,n±m g " jv"'" ^ m ."^k"' then the random matrix 
H m possesses a variance profile ll54l . Il76l . The asymptotic 
approximation for (q) is given in the following lemma. 

Lemma 1. The instantaneous random variable r{^ N (q) can be 
approximated by a deterministic quantity^ 7„ N (q) such that 



r„ (q) — 7m (q) ~~>0 as the system dimension N — >• oo. 

(q) is described by the following fixed-point equa- 



DN 



Also, 7, 
tion: 

7° N (q) 



N S 



q n q m d„ 



,d„ 



Mm. (14) 



N t-i q m d„ 



,+gnd m ,„ 7 ™(q) 



Proof: The proof is given in Appendix B. ■ 
From Lemma [T] we know that 7„ N (q) becomes more 
accurate when increasing the system dimension, and is asymp- 
totically tight for r^ N (q). For further analysis, an auxiliary 
vector 4> e R^ xl is defined with m (q) = , Vm. 

Then from Lemma Q] the fixed-point equation for </> m (q) can 
be written as 



</>m (q) 



1 



n^ra 



q n d n 



Vm. 



(15) 



>(q) 



From ( fl5l >. it is easy to see that q and <p are coupled and 
their relationship only depends on the statistical channel in- 
formation reflected in t2 m , n . Designing algorithms to compute 
q and 0(q) is of primary interest and one common approach 
is to examine the conditional convergence property of q and 
4> separately. 

The convergence property of given q is relatively easy to 
establish since it does not involve any constraint. Given any 
q satisfying the dual network power constraint, the algorithm 
to compute the corresponding 0(q) is given in Table HI] and 
is referred to as Algorithm B whose convergence property is 
given below. 



4 Note that we present the asymptotic behavior of (q) with a given 
power vector q, not with the instantaneous optimal power vector q*. 
The instantaneous optimal power vector is a function of channel and thus 
complicates standard large system analysis. Bounding techniques trying to 
investigate this issue are conducted in 1521 . In this paper, iterative method 
is used to compute the asymptotically optimal power. This comment carries 
over to the following lemmas. 
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TABLE II 

Algorithm B : Computation of given q 



TABLE III 

Algorithm C: Computation of q given <j> 



Initialize arbitrary <p[0] G 
Update 4>[t+ 1]: 

4> m [£ + l] 



with a given q. 



w m + 



E 



Vm. 



N V l + q n d„ 



Lemma 2. For a given q, starting from any initial cb\0], the 
4>[£] in Algorithm B converges to the unique solution^ of the 
fixed-point equation 1151 . 

Proof: The proof is given in Appendix B. ■ 
Now consider the convergence property of q given <p. 
Combining (fT~4T > and (Q3) yields the equivalent fixed-point 
equation for 7„ N (q) = ; — _^ mdm -™ H _ . Thus the 



E 



q n d n 



can be seen as the 



additive effect of i £ 

asymptotically equivalent interference and is regarded as ef- 
fective interference in BTI . In the sequel, we construct the 
effective dual network to draw further insight for the power 
control problem. 

Firstly, the following power control problem conditioned on 
4> is constructed by considering the weighted power constraint: 



maximize 



subject to ^■°" T q < P, 
variables : q. 



E 



q > 



Then, by defining the vector e DN = (j^- 
the nonnegative matrix E DN (q) as 



o, 



if 777 = 71 
if 772 =^ 77 



(16) 



and 



(17) 



the objective function in (fTST l can be expressed compactly as 

-, whose algebraic structure leads 



(diag(0oe°N)(^E™(q)q+w)) r 

to the following eigenvalue problem in terms of the power q* 
and weighted asymptotic SINR c;*: 



^=diag(/3oe DN ) (E 



DN 



(q*) + (l/P)w* T )^- (18) 



By comparing with (0, we can see that E DN (q) can be 
regarded as the effective cross channel interference matrix 
and the effective dual network can be characterized by the 
nonnegative matrix diag (/3 o e DN ) (E DN (q) + (l/P)wer T ). 
Note that in the finite system setting, the cross channel inter- 
ference matrix F is independent of the power. However, for 
the large system setting, E DN (q) and q are interdependent. In 
the following, we employ nonlinear Perron-Frobenius theory 
to propose a distributed algorithm to compute q* given cj>, 
which is given in Table Hill and is referred to as Algorithm C. 



5 The existence of the solution can be shown by employing the same method 
as in the proof of Theorem 6.1 in [46] and the proof of Theorem 1 in 1741 . 



. Initialize arbitrary q[0] G 

^<T T q[0] < P. 
1) Update dual network power q[£ + 1]: 

Pr, 



with a given <p such that 



Q n \$\ dm , n 



[£]dm,n4>n 



Vm. 



2) Normalize q[£+ 1]: 

q[*+l] 



NP 



<r T q[£ + 1] 



q[£+l]. 



Theorem 2. For a given <fi, starting from any initial q[0], 
the q[£] in Algorithm C converges geometrically fast to the 
optimal solution q*(0) of ( TToT l. 

Proof: The proof is given in Appendix A. ■ 
After establishing the convergence properties of computing 
<A(q) in Algorithm B and q* (<j>) in Algorithm C, both Algo- 
rithm B and Algorithm C can be combined using alternate op- 
timization to compute an asymptotically local optimal solution 
</>(q* ) and q* . The asymptotically optimal dual network power 
is used to design the asymptotically optimal beamformer 
in (fTTl i. The procedure to derive the asymptotically optimal 
primal network power requires q* and 0(q*), and is pursued 
next. 



B. Asymptotic Analysis for the Primal Network 

Similar procedure for analyzing the dual network can be 
applied to the primal network in order to examine the asymp- 
totically optimal transmit power p*. From the analysis in 
Section |nll the primal network SINR is given as 



Pmdr, 



r™(p) 



N 



Iht u* I 

I m,m m I 



E 



N 



lh T ii* I 2 
| 1A rn,m u.^ | 



(19) 



The asymptotic approximation of T™ (p) is presented in the 
following lemma. 

Lemma 3. The instantaneous random variable T™ (p) can be 
approximated by a deterministic quantity 7„ N (p) such that 
r™(p) — 7™(p) -^>0 as the system dimension N — > oo. 
Also, 7™(p) is described by the following equation: 



7™(p) = 



where <f/ m (q) = 



P 



m u "m+m . 



-0m (q) 



Vm (20) 



E 



,0 m (q)) a <(q) 



Vt77. 



Proof: The proof is given in Appendix B. ■ 
We can see from ( f20b that the effective interference 
in the primal network equals the following: 
if V tt— — 7 . In order to establish 

the effective primal network, we consider the following 
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TABLE IV 

Algorithm D: Computation of p given q and 4> 



• Initialize arbitrary p[0] £ R^^ xl with given q and cj> such that 
fW T p[0] < P. 



1) Update primal network power p[£ + 1] 

t>m[e+ 1] 



<t>? n d m , 

2) Normalize p[£ + 1]: 



"ro+iv L ~ 



^ j (l+9m<i m , m <? m ) 2 



Vm. 



AP 



w 1 p[£ + 1] 



constructed power control problem conditioned on q and <f>: 



maximize mm 



N 



(21) 



subject to jf-w 1 p < P, p > 
variables : p. 



Then, by defining the vector 

tx- . . . ILlk — i and the nonnegative matrix 



E PN as 



E 



PN 



' djK,JK<i 



0, 



if 771 = 77 

if m ^ n 



(22) 

— ii / ; , j 1 1 

the objective function in (fJTJ can be expressed compactly as 
— — — st^tt — \\ — , whose algebraic structure leads to 

(diag(/3oePN)(_LEPNp+<T)) m ' & 

the following eigenvalue problem in terms of the power p* 
and weighted asymptotic SINR £*: 



El 



diag (/3 o e PN ) (E PN + (l/P)<xw T ) (23) 



By comparing with (0, we can see that E PN can be regarded 
as the effective cross channel interference matrix and the effec- 
tive primal network can be characterized by the nonnegative 
matrix diag (/3 o e PN ) (E PN + (l/P)crw T ). Compared with 
the effective dual network, E PN is not explicitly dependent 
on p. In the following, we employ Perron-Frobenius theory 
to propose a distributed algorithm to compute p* given q and 
<j>, which is given in Table [IV] and is referred to as Algorithm 
D. 

Theorem 3. For given q and 0, starting from any initial p[0], 
the p[£] in Algorithm D converges geometrically fast to the 
optimal solution p* (q, <f>) of (T2TI) . 

Proof: The proof is given in Appendix A. ■ 
Now, by combining Algorithm B, C, and D that have 
respectively treated <f>, q and p separately, a single timescale 
algorithm is given in Table [V] and is referred to as Algorithm 
E. Even though this algorithm that computes the asymp- 
totic power is iterative, it only requires statistical channel 
information and thus the asymptotic power is updated at 
a slower timescale. Then for each instantaneous time, the 
asymptotic primal network power p* is used for the down- 
link transmission, and the asymptotic dual network power 



q* is employed to non-iteratively obtain the instantaneous 
beamforming matrix U* with local CSI as u^(q*) = 



,hi, n +w m l)~ 



— In this way, by leveraging 

!!(£„,!„, -#l»m,nhm,n+tUmI) -:l hm,mll 

the asymptotic property in the large scale system, no instan- 
taneous power update is required in the coordinated cluster 
to jointly optimize power control and beamformer. To draw 
connection with the finite system analysis, we summarize the 
results obtained by the nonlinear Perron-Frobenius theory in 
Table LYD 

Discussion of Complexity: It is important to note that even 
though Algorithm A and Algorithm E are both discrete time 
algorithms, their operating timescales as well as the imple- 
mentation complexities are vastly different (we use indices k 
and I to differentiate them). In Algorithm A, the power update 
is on the order of milliseconds to track the instantaneous 
channel effect. Thus, this algorithm requires a large amount of 
instantaneous power update to compute the optimal solution. 
In contrast, the power update in Algorithm E relies only 
on statistical channel information. Therefore, this algorithm 
operates on the order of tens of seconds or more (at the same 
timescale as the variation of the long-term channel statistics) 
and thus the implementation complexity is greatly reduced. 

V. Numerical Results 

In this section, we conduct a numerical study to support 
the analysis. We consider a three-cell cluster, i.e., J = 3, 
wherein the three base stations jointly perform power control 
and coordinated beamforming. The path loss (in dB) model 
in 11771 is assumed with 15.3 + 37.61og 10 d for distance d in 
meters and a log-normal shadowing with standard deviation 
of 8 dB is employed. The noise power spectral density is set 
to —162 dBm/Hz. The radius of each base station is set to be 
1.5 km, and a 15 dBi antenna gain is assumed. For illustration 
purpose, the coordinated cluster is subject to a total power 
constraint, i.e., w = 1, and each user possesses the same 
priority of service, i.e., (3 — 1. The total power constraint P 
is assumed to equal 10 Watt. 

Firstly, a finite system setting is considered. Each base 
station is assumed to be equipped with N = 4 antennas and 
serves K = 4 randomly located users simultaneously. For one 
channel realization, the coordinated cluster utilizes Algorithm 
A to iteratively obtain the optimal beamformer and optimal 
power. Fig. Q] demonstrates the convergence plot of the primal 
network power for different users. It is observed that for this 
channel realization, the power converges to its optimal value 
within 3 runs of iteration. This demonstrates the geometrically 
fast convergent property of Algorithm A. Extensive numerical 
evaluations show that it converges typically within 10 runs 
of iteration for different channel realizations. In Fig. |2] the 
convergence plot of the primal network weighted SINR for 
different users is shown. Since the system metric of the 
coordinated cluster is to maintain fairness across users, each 
user's optimal SINR (w = 1) would converge to the same 
value for a given channel realization, which is verified in Fig. 

12 

Next, a large system setting is considered with = 50 
and K = 40. For a given geometry, the asymptotic SINR 
of the primal network is of interest, whose convergence plot 



TABLE V 

Algorithm E: Computation of p and q for multicell downlink 



. Initialize arbitrary 0[O] £ R++ Xl , p[0] S K++ Xl , and q[0] 6 R^ xl such that iw T p[0] < P, and -^er T q[0] < P. 
1) Update dual network power q[£ + 1]: 



q m [I + 1] = \ w m + — > 

d m . m \ N ^r 1 1 + q„ 



Vm. 



2) Normalize q[^ + 1]: 

3) Update 4>[t+l\. 

4) Update (f> [t + l\: 



NP 
<t 1 q[£ + 1] 



$ m [e + i)= f m , 1 Vm. 

C(m,m 5m [« + 1J 



5) Update primal network power p [t + 1] : 



+ Vm. 



n ^ m (l+9»[f+l]dm,^ m ['+l]) 2 



* WO.1l -<J^+l]An / 1 ^ P„M<i m ,„^[Hl] 

D m ( + 1 = — (T m H > ~ ~ 

^ + l]dm,m I (l + <j m [£+l]d m:m m [£+l]) 2 0U^+l] 



6) Normalize p[f + 1]: 

w 1 p[£ + 1] 



Vm. 



TABLE VI 

Nonlinear Perron-Frobenius characterization: finite system of the dual network (second row); finite system of the primal 
network (third row); large system of the dual network (fourth row); large system of the primal network (fifth row). 



Concave Self-Mapping 


Perron Eigenvalue 


Perron Eigenvector 


diag(/3 o g(U)) (iF T (D)q(B) + w) , VU 


p (diag(/3 o g(U)) (F T (U) + (l/P)w CT T )) 


x (diag(/3 o g(U)) (F T (U) + (1/P)w<t t )) 


diag(/3og(U))(iF(U)p(U) + <T) , VU 


p (diag(/3 o g(U)) (F(U) + (l/P)<rw T )) 


x (diag(/3 o g(U)) (F(U) + (l/P)<rw T )) 


diag(/3oe DN )(iE DN (q)q + w) 


p(diag(/3oe DN ) (E DN (q) + (l/P)w CT T )) 


x (diag(/3 o e DN ) ( E ™(q) + (l/P)w CT T )) 


diag(/3oe™)(iE™p + CT ) 


p(diag(/3oe™) (E™ + (l/P) CT w T )) 


x (diag(/3 o e PN ) (e™ + (l/P)<rw T )) 



is shown in Fig. [3]by employing Algorithm E. The SINR's of 
each user are not differentiated, and uses the same line of type 
for illustration. Note that the converged value does not depend 
on the channel realization. However, it depends on the user 
geometry, namely the large scale channel effects, which means 
different user geometries would lead to different deterministic 
equivalents for the optimal SINR in the large system. Fig. |4] 
considers the use of asymptotic result. The asymptotic primal 
network power is utilized for downlink transmission, and the 
asymptotic dual network power is leveraged to non-iteratively 
determine the instantaneous beamformer. The achieved SINR's 
for different users using the determined beamformer are 
shown, along with their mean and the achieved SINR using the 
optimal beamformer obtained via Algorithm A, for one channel 
realization. It is observed that the SINR's of different users 
employing the asymptotically optimal beamformer fluctuate 
around the optimal one, with the mean close to the optimal 
SINR. Therefore, by using Algorithm E to obtain the asymptot- 
ically optimal beamformer, the max-min fairness across users 
can be achieved in the asymptotic sense. 

Finally, in Fig. [5] we consider the use of the asymptotic 



result in a finite system with N = 4 and K = 3 and demon- 
strate the comparison of the average SINR using optimal 
beamformer and the asymptotically optimal beamformer with 
respect to the variation of the total power constraint P. Herein, 
the averaging is over the user geometries, and for a given user 
geometry, different channel realizations are drawn. It can be 
seen from Fig. [5] that the performance of applying asymptotic 
result holds well for finite system setting. Accordingly, in a 
practical system with limited backhaul constraint, the asymp- 
totically optimal power and beamformer can be developed 
and leveraged to reduce the implementation complexity and 
approach the optimal performance in the asymptotic sense. 

VI. Conclusion 

In this paper, we consider a joint optimization of beamform- 
ing and power control in a coordinated multicell downlink 
and employ the max-min formulation to enforce egalitarian 
fairness across users. The network duality is interpreted via 
a nonlinear Perron-Frobenius theoretic characterization and 
utilized to design a distributed algorithm to obtain the optimal 
solution. The iterative algorithm requires instantaneous power 
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i i 1 1 1 1 1 1 1 1 

1 23456789 10 

Iteration Number 

Fig. 1. Convergence plot of the primal network power in a finite system 
setting employing Algorithm A: (TV = 4, K = 4, J = 3, P = 10 Watt). 
Different marker types represent different users. 



20 




10 20 30 40 50 

Iteration Number 



Fig. 3. Convergence plot of the primal network weighted SINR in a large 
system setting employing Algorithm E: (TV = 50, K = 40, J = 3, P = 10 
Watt). The SINR for different users are not differentiated, but use the same 
line style for illustration. 
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Iteration Number 



Each user's SINR (asymptotically optimal beamformer) 
Mean of each user's SINR (asymptotically optimal beamformer) 
- Each user's SINR (optimal beamformer) 
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Fig. 2. Convergence plot of the primal network weighted SINR in a finite 
system setting employing Algorithm A: (TV = 4, K = 4, J = 3 P = 10 
Watt). Different marker types represent different users. 



update through the limited backhaul and does not scale well 
in a large system setting. In order to design an algorithm that 
only utilizes channel statistics, we leverage random matrix 
theory to derive deterministic equivalents for the optimal SINR 
expression, and utilize the nonlinear Perron-Forbenius theory 
to establish the notion of effective network and propose a 
fast convergent algorithm. The asymptotically optimal solution 
enables a non-iterative approach to compute the instantaneous 
beamformer and thus requires no instantaneous information 
exchange across the coordinated cluster. This paper assumes 
an independent channel model and utilizes a weighted sum 
power constraint. Investigating the impact of practical issues 
on algorithm design in a large system setting, such as channel 
estimation error and per cell power constraint are interesting 
directions of future work. 



Fig. 4. The achieved primal network SINR for each individual user using the 
asymptotically optimal beamformer in a large system setting for one channel 
realization: (TV = 50, K = 40, J = 3, P = 10 Watt). The mean of the 
achieved SINR using asymptotically optimal beamformer averaged over users 
and the the achieved SINR using the optimal beamformer are illustrated for 
comparison. 



Appendix A 

Proof of Theorem\I\ The key step to the proof is to establish 
the convergence property of the dual network power q via a 
nonlinear Perron-Frobenius theory in ||36l . The relationship 
between q* and the optimal weighted SINR r* is of interest, 
and can be obtained by substituting the optimal MVDR 
beamformer: 



Nt* h' CV 2"h h 



11 m,r. 



Thus the mapping : 

can be defined by the following equation: Im\cC) — 



Vto. 
(24) 

JKxl 
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— B — Optimal beamformer 

-0- 1 Asymptotically optimal beamformer 



4 5 6 7 

Total Power (Watt) 



Fig. 5. Comparison of the average achieved primal network SINR using 
asymptotically optimal beamformer and the optimal beamformer in a finite 
system setting with respect to different values of the power constraint P. 
The averaging is performed over different geometries of users and different 
channel realizations: (N = 4, K = 3, J = 3). 



-. It can be shown using 



the same technique in [36] that I^(-) is a concave self- 
mapping of q* . Also, for the dual network, the weighted sum 
power constraint -i-er T q* = P induces a norm on M^_ Kxl 
defined by ||q*|| DN 4 (N/P)J2 m a m q* l . By applying M 
Theorem 1], starting from any initial point q[0], the fixed- 
point iteration (step 1 and step 2 of Algorithm A) converges 
geometrically fast to the optimal solution q* for the eigenvalue 
problem (f24-b . The optimal beamforming matrix U* is unique 
and can be computed by substituting the optimal dual network 
power q* into the MVDR beamformer (fTTT i for each user 
(step 3 of Algorithm A). For the primal network power p, the 
induced norm on M^ Kxl is established by the weighted power 
constraint ^w T p* = P as: ||p*||pn = (N/P) E m w m p* m - 
Therefore, by using the same line of argument for the dual 
network with the optimal beamforming matrix U*, the fixed- 
point iteration (step 4 and step 5 of Algorithm A) converges 
geometrically fast to the optimal solution p* for the eigenvalue 
problem (01 with any initial point p[0]. This completes the 
proof of Theorem [TJ 

Proof of Theorem [2} For a given <fi, the nonlinear 
eigenvalue problem in (fTol l enables us to define the map- 
ping X( 3 )(-) : R^ xl -> Kf^ 1 as: 4?(q) = 

Since the function 



(w m + Jf E 



1+x 



is strictly concave in x G the mapping Tm (q) is 

a summation of strictly concave functions in q and thus is 
a concave self-mapping in q. Then using the norm ||q|| dim 
in Appendix A and applying ||35l Theorem 1], the fixed- 
point iteration (step 1 and 2 of Algorithm C) converges 
geometrically fast to q* (</>) for the eigenvalue problem ([T8"V 

Proof of Theorem [3} For given and q, the eigen- 
value problem in ((23) enables us to define the map- 



ping K 4 )(-) 



JKxl 



JKxl 



as: l { m\f>) 



N E 



N V (l+9m<i„ 

r(4). 



It can be easily 



seen that the mapping I,„ (p) is affine, thus it is a concave 
self-mapping in p. Then using the norm ||p||pn in Appendix 
A and applying ||35l Theorem 1], the fixed-point iteration 
(step 1 and 2 of Algorithm D) converges geometrically fast 
to p* (<p, q) for the eigenvalue problem (l23l . 

Appendix B 

Useful Results from Random Matrix Theory: We reproduce 
the following theorem [48 1, [52 1, [76 1 that will be employed 
to prove Lemma Q] and Lemma [3] 

Theorem 4. (Theorem 2 in |76|) Consider an N x n random 
matrix Y = (^.j)^J™ J=1 where the entries are given by: 
Yij — ^fiXij, the Xij being independent and identically 
distributed (i.i.d.), with the following assumptions hold: 

Al: The complex random variables Xij are i.i.d. with 
E[X id ] = 0, E[Xf tj ] = 0, n\X h] \ 2 } = 1, and E[\X U \ S ] < 
oo. 



A2: There exists a real number (T n 



< oo such that: 



sup max < cr max . 

n>l l<i<N 

There exists a deterministic JVxJV matrix-valued function 
^f(z) = diag(^i(z), . . . , ipfc(z)) analytic in C — R + such 
that: 



-Ltt (YY^ — zlfr) 1 -^Tr(^(z))°^0 forzeC-R+ 
N N 

(25) 

whose elements are the unique solutions of the deterministic 
system of TV + n equations: 



-1 



-i 



for 1 < i < N 



for 1 < j < n (26) 



such that -i-Tr(\I>(z)) is the Stieltjes transform ||451 of a 
probability measure. 

Proof of Lemma Q} The technique to establish the deter- 
ministic equivalent for 7° N (q) lies in the asymptotic be- 
havior of the empirical distribution of the eigenvalue for 



(e, 



-h 

11 m,n 11 m .n 



■ W rn I 



This uplink problem 



for the equal power system has been addressed in [47], and 
the general treatment using the notion of variance profiles for 
random matrices is provided in [76 1 . Applying ||78~1 Lemma 
2.7] yields (27). 

Since the separable variance profile for the Gram ma- 
trix En^m ^'" "mA ,n is characterized by the optimal 
power q n and the large-scale channel effects e? m . n , there exist 
a deterministic equivalent for the Stieltjes transform lt45l of 
this Gram matrix. In order to invoke Theorem |4] the channel 
model needs to satisfy the two assumptions (i.e., Al and 
A2) described above. Note that the channel model in (flZb 
constitutes a special case of the channel model assumed in 
[54 1, [76], therefore the matrices considered satisfy the two 



1 1 



T& N (q)- 



Qrnd 



N 



Tr 



OVA 



-h h f 



(27) 



— |h f u*| 2 = 



N 



2 _ 



(28) 



— |h f u*| 2 

I m.rn n I 



J_ht [ V gjdnj h ht +w II h ht IV 9jd "- J h ht 4- «; I I h 



(29) 



J_ut I gj^.j u h T +?/; I I h 



J_ht | V gjrf "- 3 h h t + it; I I h IV ht +ui l| h 



(30) 



1 + q™d m , m ci 



2 SLp^jKj + w„I h„ 



necessary assumptions. Employing Theorem |4] generates the 
fixed-point equation for 7„ N (q) in (Tl4l >. 

Proof of Lemma |2} For a given q, define the following 
mapping: 1$ {<j> m ) = 1 s d — . The idea for 

proving this lemma is to use the standard interference function 
framework lED . It is straightforward to check that the posi- 
tivity and mono tonicity conditions in [21 1 hold for Xm ((f> m ). 
Also, for all e > 1, we have — , > 

- — — , which establishes the scalability con- 

dition in [21 1. Since the mapping is a standard interference 
function, the convergence result follows from lETI . thus com- 
pleting the proof of Lemma [2] 

Proof of Lemma \3][ The expression for (p) is given 
in $1% , and the optimal beamformer is the MVDR 
beamformer in (fTTl i. The asymptotic approximations for 
jf\*4n,m u m\ 2 an d j? |h}„ im u* | 2 need to be determined. The 
expression for -^|hj n m u*„| 2 can be further expanded as (28). 

Employing Theorem |U the numerator of (l28t converges 
almost surely to <^(q). In order to obtain the determin- 
istic equivalent for the denominator, the dependence of 
<f> m (q) on the noise variance w m can be made explicit, 
i.e., <j> m (q) = <f> m (q_,x)\ x=Wm . Then, by employing the 
differential of the Stieltjes transform of the Gram matrix 
E^m fe jy"" VA,n and aPPiying Theorem gj the de- 
nominator of (f28b converges almost surely to — 4 > 'm{ c L) = 
-^<^ m (q,a;)| x=Wm , which can be shown to be: 0' m (q) = 

-^m(q) 

j 1 ^ Qti dm, 77 ' 

^,77+ JV ^ (l + ,„ im ,„ +m(q ))i 

The expression for -^Ih^^u*! 2 can be further expanded 



as (29). 

Following the same line of argument, the denominator of 
d29] > converges almost surely to — <j>' n (q)- For the numerator of 
d29l , since h m ,m and the Gram matrix J2 J7 t n ^N^^^nj 
are not independent, the numerator of (|29~t is transformed into 
the equivalent form (30) by matrix inversion lemma. 

By employing the rank-1 perturbation lemma ll79l and 
Theorem the numerator of (l30t converges almost surely to 
— <fi' m (q), and the denominator of (l30l > converges almost surely 
to (1 + qmd m ,m4>m{ci)) 2 ■ Combining the aforementioned re- 
sults yields the fixed-point equation for 7™(p) in 1201 . This 
completes the proof of Lemma [3] 

Acknowledgment 

The authors would like to thank the anonymous reviewers 
and the Associate Editor for their suggestions which helped 
to improve the quality of the paper. 

References 

[I] D. Gesbert, S. Hanly, H. Huang, S. Shitz, O. Simeone, and W. Yu, 
"Multi-cell MIMO cooperative networks: a new look at interference," 
IEEE J. Sel. Areas Commun., vol. 28, no. 9, pp. 1380-1408, Dec. 2010. 

[2] J. G. Andrews, W. Choi, and R. W. Heath, "Overcoming interference in 
spatial multiplexing MIMO cellular networks," IEEE Wireless Commun., 
vol. 14, no. 6, pp. 95-104, Dec. 2007. 

[3] D. J. Love, R. W. Heath, V. K. N. Lau, D. Gesbert, B. D. Rao, and 
M. Andrews, "An overview of limited feedback in wireless commu- 
nication systems," IEEE J. Sel. Areas Commun., vol. 26, no. 8, pp. 
1341-1365, Oct. 2008. 

[4] R de Kerret and D. Gesbert, "Degrees of freedom of the network MIMO 
channel with distributed CSI," IEEE Trans. Inf. Theory, vol. 58, no. 11, 
pp. 6806-6824, Nov. 2012. 

[5] A. Lozano, R. W. Heath, and J. G. Andrews, "Fundamental limits of 
cooperation," IEEE Trans. Inf. Theory, 2013, to appear. 

[6] A. Sanderovich, O. Somekh, H. V. Poor, and S. Shamai, "Uplink macro 
diversity of limited backhaul cellular network," IEEE Trans. Inf. Theory, 
vol. 55, no. 8, pp. 3457-3478, Aug. 2009. 



12 



[7] T. Le, S. Nasseri, A. Zarrebini-Esfahani, M. Nakhai, and A. Mills, 
"Power-efficient downlink transmission in multicell networks with lim- 
ited wireless backhaul," IEEE Wireless Commun., vol. 18, no. 5, pp. 
82-88, Oct. 2011. 

[8] L. Zhou and W. Yu, "Uplink multicell processing with limited backhaul 
via successive interference cancellation," in IEEE Global Communica- 
tions Conference (Globecom), Dec. 2012, pp. 1-6. 
[9] G. Caire and S. Shamai, "On the achievable throughput of a multi- 
antenna gaussian broadcast channel," IEEE Trans. Inf. Theory, vol. 49, 
no. 7, pp. 1691-1706, Jul. 2003. 

[10] H. Weingarten, Y. Steinberg, and S. Shamai, "The capacity region of 
the gaussian multiple-input multiple-output broadcast channel," IEEE 
Trans. Inf. Theory, vol. 52, no. 9, pp. 3936-3964, Sept. 2006. 

[11] M. K. Karakayali, G. J. Foschini, and R. A. Valenzuela, "Network co- 
ordination for spectrally efficient communications in cellular systems," 
IEEE Trans. Wireless Commun., vol. 13, no. 4, pp. 56-61, 2006. 

[12] S. S. Christensen, R. Agarwal, E. Carvalho, and J. Cioffi, "Weighted 
sum-rate maximization using weighted MMSE for MIMO-BC beam- 
forming design," IEEE Trans. Wireless Commun., vol. 7, no. 12, pp. 
4792^1799, 2008. 

[13] L. Venturino, N. Prasad, and X. Wang, "Coordinated linear beamform- 
ing in downlink multi-cell wireless networks," IEEE Trans. Wireless 
Commun., vol. 9, no. 4, pp. 1451-1461, 2010. 

[14] E. Bjornson, R. Zakhour, D. Gesbert, and B. Ottersten, "Coopera- 
tive multicell precoding: rate region characterization and distributed 
strategies with instantaneous and statistical CSI," IEEE Trans. Signal 
Process., vol. 58, no. 8, pp. 4298^310, Aug. 2010. 

[15] C. W. Tan, M. Chiang, and R. Srikant, "Maximizing sum rate and 
minimizing MSE on multiuser downlink: optimality, fast algorithms and 
equivalence via max-min SINR," IEEE Trans. Signal Process., vol. 59, 
no. 12, pp. 6127-6143, Dec. 2011. 

[16] Q. Shi, M. Razaviyayn, Z. Luo, and C. He, "An iteratively weighted 
MMSE approach to distributed sum-utility maximization for a MIMO 
interfering broadcast channel," IEEE Trans. Signal Process., vol. 59, 
no. 9, pp. 4331^1340, Sept. 2011. 

[17] T. E. Bogale and L. Vandendorpe, "Weighted sum rate optimization for 
downlink multiuser MIMO coordinated base station systems: centralized 
and distributed algorithms," IEEE Trans. Signal Process., vol. 60, no. 4, 
pp. 1876-1889, Apr. 2012. 

[18] W. Utschick and J. Brehmer, "Monotonic optimization framework for 
coordinated beamforming in multicell networks," IEEE Trans. Signal 
Process., vol. 60, no. 4, pp. 1899-1909, Apr. 2012. 

[19] P. Weeraddana, M. Codreanu, M. Latva-aho, and A. Ephremides, "Mul- 
ticell MISO downlink weighted sum-rate maximization: a distributed 
approach," IEEE Trans. Signal Process., vol. 61, no. 3, pp. 556-570, 
Feb. 2013. 

[20] G. J. Foschini and Z. Miljanic, "A simple distributed autonomous power 

control algorithm and its convergence," IEEE Trans. Veh. Technol, 

vol. 42, no. 4, pp. 641-646, Nov. 1993. 
[21] R. D. Yates, "A framework for uplink power control in cellular radio 

systems," IEEE J. Set Areas Commun., vol. 13, no. 7, pp. 1341-1347, 

Sept. 1995. 

[22] E. Visotsky and U. Madhow, "Optimum beamforming using transmit 
antenna arrays," in IEEE Veh. Technol. Conf. (VTC), vol. 1, Jul. 1999, 
pp. 851-856. 

[23] F. Rashid-Farrokhi, K. J. R. Liu, and L. Tassiulas, "Transmit beamform- 
ing and power control for cellular wireless systems," IEEE J. Sel. Areas 
Commun., vol. 16, no. 8, pp. 1437-1450, Oct. 1998. 

[24] M. Bengtsson and B. Ottersten, "Optimal downlink beamforming using 
semidefinite optimization," in Annu. Allerton Conf. Commun., Contr, 
Comput, Sept. 1999, pp. 987-996. 

[25] P. Viswanath and D. N. C. Tse, "Sum capacity of the vector Gaus- 
sian broadcast channel and uplink-downlink duality," IEEE Trans. Inf. 
Theory, vol. 49, no. 8, pp. 1912-1921, Aug. 2003. 

[26] S. Vishwanath, N. Jindal, and A. Goldsmith, "Duality, achievable rates, 
and sum-rate capacity of Gaussian MIMO broadcast channels," IEEE 
Trans. Inf. Theory, vol. 49, no. 10, pp. 2658-2668, Oct. 2003. 

[27] M. Schubert and H. Boche, "Solution of the multiuser downlink beam- 
forming problem with individual SINR constraints," IEEE Trans. Veh. 
Technol, vol. 53, no. 1, pp. 18-28, Jan. 2004. 

[28] W. Yu, "Uplink-downlink duality via minimax duality," IEEE Trans. Inf. 
Theory, vol. 52, no. 2, pp. 361-374, Feb. 2006. 

[29] B. Song, R. L. Cruz, and B. D. Rao, "Network duality for multiuser 
MIMO beamforming networks and applications," IEEE Trans. Com- 
mun., vol. 55, no. 3, pp. 618-630, Mar. 2007. 

[30] H. Dahrouj and W. Yu, "Coordinated beamforming for the multicell 
multi-antenna wireless system," IEEE Trans. Wireless Commun., vol. 9, 
no. 5, pp. 1748-1759, May 2010. 



[31] A. Tolli, H. Pennanen, and P. Komulainen, "Decentralized minimum 
power multi-cell beamforming with limited backhaul signaling," IEEE 
Trans. Wireless Commun., vol. 10, no. 2, pp. 570-580, Feb. 2011. 

[32] W. Yang and G. Xu, "Optimal downlink power assignment for smart 
antenna systems," in IEEE Int. Conf. Acoust., Speech and Signal 
Process. (ICASSP), vol. 6, May 1998, pp. 3337-3340. 

[33] A. Wiesel, Y. C. Eldar, and S. Shamai, "Linear precoding via conic 
optimization for fixed MIMO receivers," IEEE Trans. Signal Process., 
vol. 54, no. 1, pp. 161-176, Jan. 2006. 

[34] Y. Huang, G. Zheng, M. Bengtsson, K. K. Wong, L. Yang, and 
B. Ottersten, "Distributed multicell beamforming with limited intercell 
coordination," IEEE Trans. Signal Process., vol. 59, no. 2, pp. 728-738, 
Feb. 2011. 

[35] U. Krause, "Concave Perron-Frobenius theory and applications," Non- 
linear Analysis, vol. 47, no. 3, pp. 1457-1466, 2001. 

[36] D. W. H. Cai, T. Q. S. Quek, and C. W. Tan, "A unified analysis of 
max-min weighted SINR for MIMO downlink system," IEEE Trans. 
Signal Process., vol. 59, no. 8, pp. 3850-3862, Feb. 2011. 

[37] D. W. H. Cai, T. Q. S. Quek, C. W. Tan, and S. H. Low, "Max-min SINR 
coordinated multipoint downlink transmission-duality and algorithms," 
IEEE Trans. Signal Process., vol. 60, no. 10, pp. 5384-5395, Oct. 2012. 

[38] C. W. Tan, S. Friedland, and S. H. Low, "Spectrum management in 
multiuser cognitive wireless networks: optimality and algorithm," IEEE 
J. Sel. Areas Commun., vol. 29, no. 2, pp. 421^130, Feb. 2011. 

[39] , "Nonnegative matrix inequalities and their application to noncon- 

vex power control optimization," SIAM Journal on Matrix Analysis and 
Applications, vol. 32, no. 3, pp. 1030-1055, 2011. 

[40] L. Zheng and C. W. Tan, "Cognitive radio network duality and algo- 
rithms for utility maximization," IEEE J. Sel. Areas Commun., vol. 31, 
no. 3, pp. 500-513, Mar. 2013. 

[41] T. L. Marzetta, "Noncooperative cellular wireless with unlimited num- 
bers of base station antennas," IEEE Trans. Wireless Commun., vol. 9, 
no. 11, pp. 3590-3600, Nov. 2010. 

[42] J. Jose, A. Ashikhmin, T. L. Marzetta, and S. Vishwanath, "Pilot 
contamination and precoding in multi-cell TDD systems," IEEE Trans. 
Wireless Commun., vol. 10, no. 8, pp. 2640-2651, Aug. 2011. 

[43] F. Rusek, D. Persson, B. K. Lau, E. G. Larsson, T. L. Marzetta, 
O. Edfors, and F. Tufvesson, "Scaling up MIMO: opportunities and 
challenges with very large arrays," IEEE Signal Process. Mag., vol. 30, 
no. 1, pp. 40-60, Jan. 2013. 

[44] F. Fernandes, A. Ashikhmin, and T. L. Marzetta, "Inter-cell interference 
in noncooperative TDD large scale antenna systems," IEEE J. Sel. Areas 
Commun., vol. 31, no. 2, pp. 192-201, Feb. 2013. 

[45] A. M. Tulino and S. Verdu, Random Matrix Theory and Wireless 
Communications. Foundations and Trends in Communications and 
Information Theory, 2004. 

[46] R. Couillet and M. Debbah, Random Matrix Methods for Wireless 
Communications. Cambridge Univ Pr, 2011. 

[47] D. N. C. Tse and S. V. Hanly, "Linear multiuser receivers: effective 
interference, effective bandwidth and user capacity," IEEE Trans. Inf. 
Theory, vol. 45, no. 2, pp. 641-657, Mar. 1999. 

[48] S. Lakshminaryana, J. Hoydis, M. Debbah, and M. Assaad, "Asymptotic 
analysis of distributed multi-cell beamforming," in IEEE Personal 
Indoor and Mobile Radio Communications (PIMRC), Sept. 2010, pp. 
2105-2110. 

[49] H. Huh, A. M. Tulino, and G. Caire, "Network MIMO with linear zero- 
forcing beamforming: large system analysis, impact of channel esti- 
mation and reduced-complexity scheduling," IEEE Trans. Inf. Theory, 
vol. 58, no. 5, pp. 2911-2934, May 2012. 

[50] S. Wagner, R. Couillet, M. Debbah, and D. Slock, "Large system 
analysis of linear precoding in correlated MISO broadcast channels 
under limited feedback," IEEE Trans. Inf. Theory, vol. 58, no. 7, pp. 
4509^4537, Jul. 2012. 

[51] C. B. Peel, B. M. Hochwald, and A. L. Swindlehurst, "A 
vector-perturbation technique for near-capacity multiantenna multiuser 
communication-part I: channel inversion and regularization," IEEE 
Trans. Commun., vol. 53, no. 1, pp. 195-202, Jan. 2005. 

[52] R. Zakhour and S. V. Hanly, "Base station cooperation on the downlink: 
large system analysis," IEEE Trans. Inf. Theory, vol. 58, no. 4, pp. 2079- 
2106, Apr. 2012. 

[53] , "Min-max power allocation in cellular networks with coordinated 

beamforming," IEEE J. Sel. Areas Commun., vol. 31, no. 2, pp. 287-302, 
Feb. 2013. 

[54] W. Hachem, P. Loubaton, and J. Najim, "Deterministic equivalents 
for certain functionals of large random matrices," Ann. Appl. Probab., 
vol. 17, no. 3, pp. 875-930, 2007. 

[55] S. P. Boyd and L. Vandenberghe, Convex Optimization. Cambridge 
Univ Pr, 2004. 



13 



[56] Y. Huang, C. W. Tan, and B. D. Rao, "Asymptotic analysis of max- 
min weighted SINR in multiuser downlink," in Conf. Inf. Sciences Syst. 
(CISS), Mar. 2012, pp. 1-5. 

[57] A. Berman and R. J. Plemmons, Nonnegative Matrices in the Mathe- 
matical Sciences. New York: Academic, 1979. 

[58] S. Boyd, A. Ghosh, B. Prabhakar, and D. Shah, "Randomized gossip 
algorithms," IEEE Trans. Inf. Theory, vol. 52, no. 6, pp. 2508-2530, 
Jun. 2006. 

[59] Y. Huang, "Feedback, power control, and beamforming: methods for 
situational aware wireless networks," Ph.D. dissertation, University of 
California, San Diego, 2012. 

[60] Y. Huang and B. D. Rao, "Awareness of channel statistics for slow cyclic 
prefix adaptation in an OFDMA system," IEEE Wireless Coimnun. Lett. , 
vol. 1, no. 4, pp. 332-335, Aug. 2012. 

[61] , "Performance analysis of heterogeneous feedback design in an 

OFDMA downlink with partial and imperfect feedback," IEEE Trans. 
Signal Process., vol. 61, no. 4, pp. 1033-1046, Feb. 2013. 

[62] S. Yarkan and H. Arslan, "Exploiting location awareness toward im- 
proved wireless system design in cognitive radio," IEEE Commun. Mag. , 
vol. 46, no. 1, pp. 128-136, Jan. 2008. 

[63] S. Ghosh, B. D. Rao, and J. R. Zeidler, "Outage-efficient strategies for 
multiuser MIMO networks with channel distribution information," IEEE 
Trans. Signal Process., vol. 58, no. 12, pp. 6312-6324, Dec. 2010. 

[64] Y. Huang, C. W. Tan, and B. D. Rao, "Outage balancing in multiuser 
MISO networks: network duality and algorithms," in IEEE Global 
Communications Conference (Globecom), Dec. 2012, pp. 1-6. 

[65] W. L. Li, Y. Shen, Y. J. Zhang, and M. Z. Win, "Robust power allocation 
for energy-efficient location-aware networks," IEEE/ACM Trans. Netw., 
2013, to appear. 

[66] S. R. Cho and W. Choi, "Energy-efficient repulsive cell activation for 
heterogeneous cellular networks," IEEE J. Set Areas Commun., vol. 31, 
no. 5, May 2013. 

[67] Y. Huang and B. D. Rao, "An analytical framework for heterogeneous 
partial feedback design in heterogeneous multicell OFDMA networks," 
IEEE Trans. Signal Process., vol. 61, no. 3, pp. 753-769, Feb. 2013. 

[68] M. Hong, R. Sun, H. Baligh, and Z. Luo, "Joint base station clustering 
and beamformer design for partial coordinated transmission in hetero- 
geneous networks," IEEE J. Sel. Areas Commun., vol. 31, no. 2, pp. 
226-240, Feb. 2013. 

[69] H. Hu, J. Zhang, X. Zheng, Y. Yang, and P. Wu, "Self-configuration 
and self-optimization for LTE networks," IEEE Commun. Mag., vol. 48, 
no. 2, pp. 94-100, Feb. 2010. 

[70] O. G. Aliu, A. Imran, M. A. Imran, and B. Evans, "A survey of self 
organisation in future cellular networks," IEEE Communications Surveys 
& Tutorials, vol. 15, no. 1, pp. 336-361, Feb. 2013. 

[71] J. Dumont, W. Hachem, S. Lasaulce, P. Loubaton, and J. Najim, "On 
the capacity achieving covariance matrix for Rician MIMO channels: 
an asymptotic approach," IEEE Trans. Inf. Theory, vol. 56, no. 3, pp. 
1048-1069, Mar. 2010. 

[72] G. Taricco and E. Riegler, "On the ergodic capacity of correlated Rician 
fading MIMO channels with interference," IEEE Trans. Inf. Theory, 
vol. 57, no. 7, pp. 4123^1137, Jul. 2011. 

[73] C. Artigue and P. Loubaton, "On the precoder design of flat fading 
MIMO systems equipped with MMSE receivers: a large-system ap- 
proach," IEEE Trans. Inf. Theory, vol. 57, no. 7, pp. 4138—4155, Jul. 
2011. 

[74] R. Couillet, M. Debbah, and J. W. Silverstein, "A deterministic equiv- 
alent for the analysis of correlated MIMO multiple access channels," 
IEEE Trans. Inf. Theory, vol. 57, no. 6, pp. 3493-3514, Jun. 2011. 




[75] J. Zhang, C. K. Wen, S. Jin, X. Gao, and K. K. Wong, "On capacity 
of large-scale MIMO multiple access channels with distributed sets of 
correlated antennas," IEEE J. Sel. Areas Commun., vol. 31, no. 2, pp. 
133-148, Feb. 2013. 

[76] A. Kammoun, M. Kharouf, W. Hachem, and J. Najim, "A central 
limit theorem for the SINR at the LMMSE estimator output for large- 
dimensional signals," IEEE Trans. Inf. Theory, vol. 55, no. 11, pp. 5048- 
5063, Nov. 2009. 

[77] 3GPP, "TR 36.814 V9.0.0 - Further Advancements for E-UTRA Phys- 
ical Layer Aspects (Release 9)," 3GPP, Tech. Rep., 2010. 

[78] Z. D. Bai and J. W. Silverstein, "No eigenvalue outside the support of 
the limiting spectral distribution of large dimensional sample covariance 
matrices," Ann. Probab., vol. 26, no. 1, pp. 316-345, 1998. 

[79] J. W. Silverstein and Z. D. Bai, "On the empirical distribution of 
eigenvalues of a class of large dimensional random matrices," Journal 
of Multivariate Analysis, vol. 54, no. 2, pp. 175-192, 1995. 

Yichao Huang (S' 10-M' 12) received the B.Eng. de- 
gree in information engineering with highest honors 
from the Southeast University, Nanjing, China, in 
2008, and the M.S. and Ph.D. degrees in electrical 
engineering from the University of California, San 
Diego, La Jolla, in 2010 and 2012, respectively. He 
then joined Qualcomm, Corporate R&D, San Diego, 
CA, as a senior system engineer. 

He interned with Qualcomm, Corporate R&D, 
San Diego, CA, during summers 2011 and 2012. 
He was with California Institute for Telecommuni- 
cations and Information Technology (Calit2), San Diego, CA, during summer 
2010. He was a visiting student at the Princeton University, Princeton, NJ, 
during spring 2012. His research interests include communication theory, opti- 
mization theory, wireless networks, and signal processing for communication 
systems. 

Dr. Huang received the President Fellowship both in 2005 and 2006, and 
the Best Thesis Award in 2008 from the Southeast University. He received the 
Microsoft Young Fellow Award in 2007 from Microsoft Research Asia and 
the ECE Department Fellowship from the University of California, San Diego 
in 2008 and was a finalist of Qualcomm Innovation Fellowship in 2010. 



Chee Wei Tan (M'08-SM'12) received the M.A. 
and Ph.D. degree in electrical engineering from 
Princeton University, Princeton, NJ, in 2006 and 
2008, respectively. 

He is an Assistant Professor at City University 
of Hong Kong. Previously, he was a Postdoctoral 
Scholar at the California Institute of Technology 
(Caltech), Pasadena, CA. He was a Visiting Faculty 
at Qualcomm R&D, San Diego, CA, in 2011. His 
research interests are in wireless and broadband 
communications, complex networks, signal process- 
ing and nonlinear optimization. 

Dr. Tan currently serves as an Editor for the IEEE TRANSACTIONS ON 
COMMUNICATIONS. He was the recipient of the 2008 Princeton University 
Wu Prize for Excellence and the 2011 IEEE Communications Society AP 
Outstanding Young Researcher Award. He was a selected participant at 
the U.S. National Academy of Engineering China-America Frontiers of 
Engineering Symposium in 2013. 




< 



14 




Bhaskar D. Rao (S'80-M'83-SM'91-F'0O) re- 
ceived the B.Tech. degree in electronics and elec- 
trical communication engineering from the Indian 
Institute of Technology, Kharagpur, India, in 1979 
and the M.S. and Ph.D. degrees in electrical engi- 
neering from the University of Southern California, 



Los Angeles, in 1981 and 1983, respectively. 
M Since 1983, he has been with the University of 

Mm California San Diego, La Jolla, where he is cur- 
iQ^Rl rently a Professor with the Electrical and Computer 

Engineering Department and holder of the Ericsson 
endowed chair in wireless access networks. He was the Director of the Center 
for Wireless Communications from 2008-2011. His interests are in the areas 
of digital signal processing, estimation theory, and optimization theory, with 
applications to digital communications, speech signal processing, and human- 
computer interactions. 

Dr. Rao was elected IEEE Fellow in 2000 for his contributions in high 



resolution spectral estimation. His research group has received several paper 
awards. His paper received the best paper award at the 2000 Speech Coding 
Workshop and his students have received student paper awards at the 2005 and 
2006 International Conference on Acoustics, Speech and Signal Processing 
as well as the best student paper award at the 2006 Neural Information 
Processing Systems (NIPS). A paper he co-authored with B. Song and R. 
Cruz received the 2008 Stephen O. Rice Prize Paper Award in the Field 
of Communications Systems. Another paper, he co-authored with David 
Wipf, received the 2012 signal processing society best paper award. He 
has been a member of the Statistical Signal and Array Processing technical 
committee, the Signal Processing Theory and Methods Technical Committee, 
the Communications Technical Committee of the IEEE Signal Processing 
Society. He has also served on the editorial board of the EURASIP Signal 
Processing Journal. He is a member of the signal processing society technical 
committee on machine learning for signal processing. 



